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1 Introduction 

Duality rotations for massless spin 2 fields have recently been shown to be symmetries of 
the action in the context of a Hamiltonian formulation involving suitable potentials that 
arise when solving the constraints HJ. This symmetry can be extended to higher spin 
fields [2] and also to the case of spin 2 fields propagating on an (A)dS background [[31 IH, 
but does not survive gravitational self-interactions [5]. 

The spin 2 result generalizes the so-called double potential formalism for spin 1 fields 
ll6l 13, which has been extended so as to include couplings to dynamical dyons by us- 
ing Dirac strings [HI HI. In the original double potential formalism, Gauss's constraint is 
solved in terms of new transverse vector potentials for the electric field so that electro- 
magnetism is effectively formulated on a reduced phase space with all gauge invariance 
eliminated. Alternatively, one may choose [fTOll to double the gauge redundancy of stan- 
dard electromagnetism by using a description with independent vector and longitudinal 
potentials for the magnetic and electric fields and 2 scalar potentials that appear as La- 
grange multipliers for the electric and magnetic Gauss constraints. In this framework, the 
string-singularity of the solution describing a static dyon is resolved into a Coulomb-like 
solution. Furthermore, magnetic charge no longer appears as a topological conservation 
law but as a surface charge on a par with electric charge. 

The aim of the present work is to apply the same strategy to the spin 2 case. Dou- 
bling the gauge invariance by keeping all degrees of freedom of symmetric tensors now 
leads to a second copy of linearized lapse and shifts as Lagrange multipliers for the new 
magnetic constraints. As a consequence, the string singularity of the gravitational dyon, 
the linearized Taub-NUT solution is resolved and becomes Coulomb-like exactly as the 
purely electric linearized Schwarzschild solution. Furthermore, as required by manifest 
duality, magnetic mass, momentum and Lorentz charges also appear as surface integrals. 

Our work thus presents a manifestly duality invariant alternative to [1 IJ where the 
coupling of spin 2 fields to conserved electric and magnetic sources has been achieved in 
a manifestly Poincare invariant way through the introduction of Dirac strings. 

Recent and not so recent related work includes for instance [fT2l[T3l[T4l[T5l[T6l[T7l[T8l 
[I9l|20l|2ll|22ll23ll2l and references therein. 

In section |2l we briefly recall the two ingredients needed for our formulation: the 
Hamiltonian description of spin 2 fields propagating on Minkowski spacetime and the 
decomposition of symmetric tensors into their irreducible components, giving rise to the 
reduced phase space description of linearized gravity. 

Our analysis starts in section |3] with a degree of freedom count that shows that the 
phase space of duality invariant spin 2 fields with doubled gauge invariance can be taken 
to consist of 2 symmetric tensors, 2 vectors and 2 scalars in 3 dimensions. We then 
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define the metric, extrinsic curvature and their duals in terms of the phase space variables 
and propose our duality invariant action principle with enhanced gauge invariance. We 
proceed by identifying the canonically conjugate pairs and discuss the gauge structure, 
Hamiltonian and duality generators of the theory. In the absence of sources, we then show 
how the generators for global Poincare transformations of Pauli-Fierz theory, reviewed in 
detail in appendix |Al can be extended to the duality invariant theory. 

The coupling to external electric and magnetic sources is discussed in section HI The 
equations of motion are first solved in the simplest case of a point-particle dyon sitting 
at the origin. They are Coulomb-like without string singularities. By identifying the 
Riemann tensor in terms of the canonical variables and computing it for this case, we 
show in appendix IE that this solution indeed describes the linearized Taub-NUT solution. 

In section |5] we discuss the surface charges of the theory and show that they include 
electric and magnetic energy-momentum and angular momentum. Because of the non- 
locality of the Poisson structure, we proceed indirectly and show that the expressions 
obtained by generalizing the surface charges of Pauli-Fierz theory in a duality invariant 
way fulfill the standard properties. Finally, we investigate how the surface charges trans- 
form under a global Poincare transformation of the sources. 

2 Preliminaries 

2.1 Canonical formulation of Pauli-Fierz theory 

The Hamiltonian formulation of general relativity linearized around flat spacetime is 



SpF[hmn. vr™", n^, n] = Jdt[J d\ (vr"^"/^^, - rTU^ - nU) - Hpp 



(2.1) 



with 



1 1 

TT VTmr, - -TT + -O n Orflmn- 



2 4 

^d^h^"d'hrn + U^hd^h^^ - ^-d^hd^h) , (2.2) 



and 



nm = -2d^n^^, 7^x = A/i-9™9"/i™„. (2.3) 



Here, indices are lowered and raised with the flat space metric 6mn and its inverse, h = 
h^'rn, = 7r™m ^ud A = dmd'^ is the Laplacian in flat space. The linearized 4 metric is 
reconstructed using hoo = —2n and hoi = rii. 
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2.2 Decomposition of symmetric rank two tensors 

Symmetric rank two tensors (prnn decompose as ||25ll26ll 

0mn = 0mn + + 0mn, (2-4) 

0mn = dmi'n + dni'm, (2.5) 
4>mn = - (5mnA - (9m9„) (2.6) 

Here is the trans verse-traceless part, containing two independent components. The 
tensor contains the trace of the transverse part of (prnn and only one independent 
component. The last three components are the longitudinal part contained in In 
terms of the original tensor (pmn the potentials for the longitudinal part and the trace are 
given by 

V = A-i (0 - A-i9™9"0^„) , (2.8) 
while the transverse-traceless part is then defined as the remainder, 



^mn ~ 4'mn ~ 4'mn ~ 'Pmn- 



(2.9) 



This implies 

A20^^ = A20^„ - Aa^^Vfcn - A9„9Vfcm 

- ^A(5^„A - dmdn)<P + ^((5m„A + dmdn)d^d^<Pkh (2.10) 
j ci3a;0™"A20^^ = j d'x (A0"^"A0^„ + 29„0™"A9Vfc„ - ^{A^f 

+ 0M-A<p + l0^0.r^0.0^<P^^). (2.11) 

Alternatively, one can introduce the local operator V^^ 

i^^^^)mn = I bmp,d^i^^\ ' ^n^.^^^ + e„p,9P(A0'?„ - Omdr^'l] , (2-12) 

which projects out the longitudinal and trace parts and onto a transverse-traceless tensor, 

(P^^0) =pTT/^TTN ^r^TT^^^TT^ 

\ ^ I mn \^ / mn \ ^ I ran ^ 

In addition, 

(P^^(P^^0))_ = -A^^Z- (2.14) 
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As a consequence, the transverse-traceless tensor can be written as V^'^ acting on a 
suitable potential 



= , V^^I = -A-' {^^^<^)mn ■ (2-15) 

The operator is related to the way the Hamiltonian constraint 7i = is solved by 
expressing the metric hmn in terms of superpotentials in [HI. When acting on a transverse- 
traceless tensor, the two terms of (|2.12l) involving 9^0'^'' can be dropped. In this case, 'P-'"^ 
is related to the generalized curl [l2l[5l, 

W)^n = \{^m,,d^'^\ + 6.,,a^0^^), (2.16) 

I^TT^TT\ =^(OfT\ . (2.17) 

\ ^ I ran \ ^ I ran ^ ' 

A second operator that projects out the longitudinal and trace parts and onto a transverse- 
traceless tensor is Q^^, 

= er.p,en..aPa'-A0''^ - ^((5„„A - 9^9„)(A0 - a^S^p,). (2.18) 

In this case, 

= AV^^, (2.19) 



so that the transverse-traceless tensor can be written as acting on another po- 
tential 



'^^-(S^V)„„. X-JJ„ = A-(S^^*)„„. (2.20) 



In turn this operator is related to the way the constraints Ti^ = are solved by expressing 
the momenta vr™" in terms of superpotentials in When acting on a transverse-traceless 
tensor, the last term can again be dropped and it is related to the square of the generalized 
curl, 

(2""'^"").. = A {0{0f-))_ = -A^G. (2.21) 

The elements of the decomposition are orthogonal under integration if boundary terms 
can be neglected, 

d X(p (fimn = d X[(p + (p + (p . (2.22) 

and the operators V^^, Q'^^, O are self-ajoint, e.g., 

j d'x (P^^0) '"" = I d'xr" {V^\)^r. ■ (2.23) 
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2.3 Reduced phase space for linearized gravity 

For completeness, let us briefly recall f25\ the reduced phase space associated with the 
Pauli-Fierz action. Because of the orthogonality of the decomposition, the canonically 
conjugate pairs can be directly read off from the kinetic term and are given by 

{hZ{x),7r'My)), Knix),<{y)), (/^L(^),4^(y))- (2.24) 

The first class constraints Hm = = H are equivalent to tt^' = = They can 

be gauge fixed through the conditions = = nff. The reduced theory only de- 
pends on 2 degrees of freedom (per spacetime point), the transverse-traceless components 
{h'^{x), T^Triy)) ^iid the reduced Hamiltonian simplifies to 

H^ = j d'x (nJp^nZ + IdrhZd'h^^) . (2.25) 

3 Action and symmetries 
3.1 Degree of freedom count 

In order to be able to couple to sources of both electric and magnetic type in a duality 
invariant way, we want to keep all components and double the gauge invariance of the 
theory. With 2 degrees of freedom, ^ dof = 2, and 8 first class constraints, # fee = 8, we 
thus need 10 canonical pairs, # ep = 10, according to the degree of freedom count ll27l 

2* (#ep) = 2* (#dof) + 2* (#fee). (3.1) 

This can be done by taking 2 symmetric tensors, 2 vectors and 2 scalars as fundamental 
canonical variables. 



A _ (Tj-a Aa (^a\ 



(3.2) 
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3.2 Change of variables and duality rotations 

For a = 1, 2, consider = /?.^„) and vr™" = (vr^", vr™-") and the definitions 

+ drX^-^e„,„d^drH-'i' + A^) + ^(5„„A - (3.3) 
= A-i - d^md'^H:^^ - dnd'-H^^ - ^(5^„A - d^d^W 

= -^H^^n- (3.4) 

The relations for hmn[H^, A^, C^] and 7r'""[i7^] are the local change of coordinates from 
the standard canonical variables of linearized gravity to the new variables. They are 
invertible and, as usual, the inverse is not local. The relations for h"^^ = h^^, n^"- = 
^mn ggj^e denote convenient combinations of the new variables in terms of which 
expressions below will simplify. As indicated by the notation, the infinitesimal duality 
rotations among the fundamental variables are 

Here, eah is skew-symmetric with ei2 = 1 and indices are lowered and raised with 5ab 
and its inverse. Duality invariance will be manifest if all the internal indices a are con- 
tracted with the invariant tensors 5ab, ^ab- Since vr™" are linear combinations of the 
fundamental variables, they are rotated in exactly the same way. We can thus consider 
Knn = ^mn' '^T^ = ^d" ^^c dual Spatial metric and the dual extrinsic curvature in the 
linearized theory. 



3.3 Action principle and locality 

The duality invariant local action principe that we propose is of the form 

Sg[z^,u''] = j d^iaAW^-u'^iM)- j dtH[z], (3.6) 

where m° denote the 8 Lagrange multipliers and 7q, the constraints. 

Let us stress here that we use the assumption that the flat space Laplacian A is invert- 
ible in order to show equivalence with the usual Hamiltonian or covariant formulation 
of Pauli-Fierz theory and also to disentangle the canonical structure. The action prin- 
ciple (|3.6I) itself and the associated equations of motion will be local both in space and 
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in time independently of this assumption. The theory itself is not local as a Hamiltonian 
gauge theory (see e.g. If28l . chapter 12) because the Poisson brackets among the canonical 
variables will not be local. 



3.4 Canonical structure 

The explicit expression for the kinetic term is 

+ ^(5™" A - 9™a")AC'^) . (3.7) 
The canonically conjugate pairs are identified by writing the integrated kinetic term as 

- 2Ad^Hl^-Al - ^A(A/7| - dpd,Hl^'')C^ + ^ A(Ai7i - dpd,H'^'^)C''') . (3.8) 

This means that the usual canonical pairs of linearized gravity can be choosen in terms of 
the new variables as 



yZ^ix), -2A {OH'^^f (y)), [C\x), -^A(Ai7| - d,d,H'/'^) (y)] 

[Al{x),2AdrHl^-{y)), (3.9) 

The 4 additional canonical pairs are 

(^Alix), -2AdrHl^^{y)), {c\x), \a{AH't - d^dXr^'M)- (3-10) 



In particular, it follows that 



{hZi^), n"^\y)} = e'^'liS'jl + S'jl)6^'\x, y). 



(3.11) 



3.5 Gauge structure 

The constraints 7q, = {Ham, 'Hai_) are choosen as 



'Ham = —'^^abd^T^mn = SCafc A(9"if^„, (3.12) 

.2, 



Hal. = Aha-dmdnhr = ^'Ca. (3.13) 

They are first class and abelian 

{7a,7/3} = 0. (3.14) 
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The constraints Him,Hi± are those of the standard Hamiltonian formulation of Pauli- 
Fierz theory expressed in terms of the new variables. The new constraints 7^ = are 

n2m = o = n2±. (3.15) 

They are equivalent to d^H^^^ = = and are gauge fixed through the conditions = 
= H^^. This does not affect vr^^', while is changed by a gauge transformation. The 
partially gauge fixed theory corresponds to the usual Pauli-Fierz theory in Hamiltonian 
form as described in section llHl 

More precisely, the observables of a Hamiltonian field theory with constraints are 
defined as equivalence classes of functionals that have weakly vanishing Dirac brackets 
with the constraints and where two functionals are considered as equivalent if they agree 
on the surface defined by the constraints (see e.g. [|28l ). The new constraints together with 
the gauge fixing conditions form second class constraints. The Dirac bracket algebra 
of observables of this (partially) gauge fixed formulation is isomorphic to the Poisson 
bracket algebra of observables of the extended formulation on the one hand, and to the 
Poisson bracket algebra of observables of Pauli-Fierz theory on the other hand. 

In the same way, the original constraints Him = = Hi^ are equivalent to -f^^m — 
= and are gauge fixed through = = H'^^, leading to the completely reduced 
theory in terms of the 2 transverse-traceless physical degrees of freedom. 

If = (^"'",^""'") collectively denote the gauge parameters, the gauge symmetries 
are canonically generated by the smeared constraints, 

S,z^ = {z^,r[e]}, Tie] = J ci^x7„e", (3.16) 

so that 

6,H^^ = -A-'e-\6mnA-dmdn)^t, = C 5^^" = 0, (3.17) 

which implies in particular 

S^hl^ = + dnCm, 4vr»„ = e-\6"^^A - 0^8^)^^. (3.18) 

Note that a way to get local gauge transformations for the fundamental variables is to 
multiply the constraints by A, which is allowed when the flat space Laplacian is invertible. 
This amounts to introducing suitable potentials for the gauge parameters and Lagrange 
multipliers. 
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3.6 Duality generator 

The canonical generator for the infinitesimal duality rotations (13.51) is 

D = jd'x(^- {V^^H-)mnHr + 2^drH:^Al^ 

~ _ y d?xV'^^{H'^)mnHr- (3.19) 
The duality generator is only weakly gauge invariant, 

{Ham-iD} = tabH^n {HaL, D} = eabH^. (3.20) 

On the constraint surface, it coincides with the generator found in pj up to normalisation, 
where it has been cast in the form of a Chem-Simons term. 

3.7 Hamiltonian 

In terms of the new variables (I3.3I) - (|3.4I) . the Pauli-Fierz Hamiltonian reads 

HpF = J d^x (if"""A2i7jJ„ - 2/\d'HlJ,H^ 



r2sn 



- d'd'H^.AH^ - Ud'd'Hlf + IaC^A^C") , (3.21) 
2 8 / 

where one can use (I2.10|) to expand the first term as a local functional of H'^^. 

The local Hamiltonian H = J (fx h of the manifestly duality invariant action princi- 
ple (Il6l) is 



d'x (^AH^^AHr - \^H-^H, + i AC^^A^a) • (3.22) 

It is equivalent to the Pauli-Fierz Hamiltonian since it reduces to the latter when the addi- 
tional constraints d'^'Hl^ = = hold. Note that the terms proportional to d' H^^ and 
may be dropped since they vanish on the constraint surface, H ^ J d^x if"™" A^if^„. 

The Hamiltonian is gauge invariant on the constraint surface, 

{H,Tm = J d^xUld^il (3.23) 

In order for the action (13.61) to be gauge invariant, it follows from (13.231) that the 
Lagrange multipliers need to transform as 

= ^'^^ - d'^C^, d^u"^ = C^. (3.24) 
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3.8 Poincare generators 

Consider now a symmetry generator of Pauli-Fierz theory. It is defined by an observable 
K[h^, TT^] whose representative is weakly conserved in time, 

^^K + {K,Hpf}^0. (3.25) 

Since the new Hamiltonian differs from the Pauli-Fierz one by terms proportional to the 
new constraints = given explicitly in (13.151) . we have 

H = HpF + j d^xj^k^. (3.26) 

Furthermore, since K, when expressed in terms of the new variables, does not depend on 
Hi, H)p, so that {K, J (Px'y^k'^} in the extended theory, it follows that K is also 
weakly conserved and thus a symmetry generator of the extended theory, 

^^K + {K,H}^0. (3.27) 

Consider then the Poincare generators Qg{^, o.) of Pauli-Fierz theory as described in 
Appendix |A1 When expressed in terms of the new variables, they are representatives for 
the Poincare generators of the extended theory. Indeed, we just have shown that they are 
symmetry generators, while we have argued in Section |33] that their Poisson algebra is 
isomorphic when restricted to their respective constraint surfaces. 

Since symmetry generators form a Lie algebra with respect to the Poisson bracket, 
{Qduj, a), D} is a symmetry generator. In much the same way as for the Hamiltonian, 
we now want to show that one can find representatives for the Poincare generators that 
are duality invariant, 

{Q^{u,a),D} = 0, (3.28) 
by adding terms proportional to the new constraints. 

The first step in the proof consists in showing that the reduced phase space generators, 
i.e., the generators Qai^^, o.) for which all variables except for the physical H^q. have been 
set to zero, are duality invariant. All other contributions to Qg{^, a) are then shown to 
be proportional to the constraints of Pauli-Fierz theory. Both these steps follow from 
straightforward but slightly tedious computations. For the generators of rotations and 
boosts for instance the computation is more involved because the explicit dependence 
has to be taken into account when performing integrations by parts. 

In terms of the new variables, the terms proportional to the constraints are bilinear 
in {h^,A'^), {tt^,A'^), {h^,C^) and {n^,C^). The duality invariant generators (5g(u;,a) 
are then obtained by adding the same terms with the substitution h'^, —A^, 

vr^ —71^ and — > C^, while keeping unchanged the terms involving only the physical 
variables H^j,. 
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As a consequence, the duality invariant Poincare transformations of /i^,7r^ are un- 
changed on the extended constraint surface. They are given by (|A.31I) - (|A.32I) where 

= -J^yx"" + a° and = -to'^x" + a\ Because of (l3^ . those of h"^, -n^ are 
obtained, on the contraint surface, by applying a duality rotation to the right hand- sides 
of (lA3T]) - (lA32l) . 

An open question that we plan to address elsewhere is the construction of the canoni- 
cal generators for global Poincare transformations in the presence of both types of sources 
that will be introduced in the next section. 

4 Coupling to conserved electric and magnetic sources 

4.1 Interacting variational principle 

We define 

^Om — — KnOi ^00 ~ ~2n", (4.1) 

and consider the action 

St[z^, u^- T"^1 = j^Sg + S', (4.2) 
IdttG 

with Sg given in (|3.6I) and the gauge invariant interaction term 

sJ = j d^x h^rr, d,Tr = 0, (4.3) 

where T^^ = (T^*^, 6'^'') are external, conserved electric and magnetic energy-momentum 
tensors. 

4.2 Linearized Taub-NUT solution 

We start by considering the sources corresponding to a point-particle gravitational dyon 
with electric mass M and magnetic mass N at rest at the origin of the coordinate system, 
for which 

rrix) = 5^^5lMJ^^\x'), Ma = (M, N). (4.4) 
In this case, only the constraints (|3.13l) are affected by the interaction and become 

HaL = -l67rGMJ^^\x). (4.5) 



They are solved by 



AC = GM^i-), 
r 



(4.6) 
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where r = ^x^Xi. It is then straightforward to check that all equations of motions are 
solved by 

= GM%2r), = G'M"(--), = n"'™ = H"^^ = 0, 

r 

/i;;„ = GM'^(<5_ + ^), vrr = 0. (4.7) 

The usual Schwarzschild form is obtained after a gauge transformation with parameter 
C'" = GM''{-1^), = 0. The solution then reads 

= GM\2r), n" = GM"(-i), A"^ = GiVrf--— ), n""^ = H^^ = 0, 

r 2 r 

= GM'^(^), nr = 0. (4.8) 

By computing the Riemann tensor in terms of the canonical variables in Appendix E 
we show that this solution describes the linearized Taub-NUT solution. It resolves the 
string singularity of the linearized Taub-NUT solution in the standard Pauli-Fierz formu- 
lation. In spherical coordinates, the latter can for instance be described by 

2GM 

hrr = = hoo, ho^ = -2N{1 - COS 9), (4.9) 

and all other components vanishing, with a string-singularity along the negative 2-axis. 



5 Surface charges 

5.1 Regge-Teitelboim revisited 

Because the theory is not local as a Hamiltonian gauge theory, the analysis of surface 
charges cannot directly be performed as in [|29l[30ll . We thus revert to the original Hamil- 
tonian method of [|3n[32ll and adapt it to the present situation of exact solutions, where 
there is no need to discuss fall-off conditions. 

Let Ch = claz"^ — h — 'jau"', with h a first class Hamiltonian density and 7q, first 
class constraints and define 0* = (z^,u"). Even though it is not so for our theory, let 
us first run through the arguments in the case where one has Darboux coordinates for the 
symplectic structure, i.e., when aAB = — is the constant symplectic matrix. We 
furthermore suppose that we are in a source-free region of spacetime. In this case one can 
show that 

4^^^ + = -ao(7a^") - d.sl, (5.1) 

where = sl[z, u] vanishes when the Hamiltonian equations of motion, including con- 
straints, are satisfied, s* ^ 0. This identity merely expresses the general fact that the 
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Noether current associated to a gauge symmetry can be taken to vanish when the equa- 
tions of motions hold (see e.g. Il28l . chapter 3), ~ 0, and that the integrand of the 
generator is given by (minus) the constraints contracted with the gauge parameters in 
the Hamiltonian formalism, s° = — 7q,£:". An explicit expression for in terms of the 
structure functions can for instance be found in Appendix D of \30i . Using integrations 
by parts, one can write the variations of the constraints under a change of the canonical 
coordinates z"^ as an Euler-Lagrange derivative, up to a total derivative, 

5,(7«£°) = 5;2^%P -9,A;:. (5.2) 

where fc* = kl[5z, z] depends linearily on and its spatial derivatives. Taking the time 
derivative of (|5.2I) and using a variation 6^ of (15.11) to eliminate dodzi'ja^"), one finds 

d.{dokl - 5,4) = 5o(fc^%|^) + S,[6,z^'-^ + S.u-'-^). (5.3) 

One now takes e° to satisfy 6i;^zf = = S^^u". Note that in the case of Darboux 
coordinates, this also implies that ^"a' = 0. If furthermore z^, is a solution of the 
Hamiltonian equations of motion, the RHS of (15.31 ) also vanishes. By using a contracting 
homotopy with respect to 50' and their spatial derivatives, one deduces that 

dokl[6z,z,] = i5^sl)[z,]-d,kl^\ (5.4) 

where A;!*/' = k£^\6(f), (pg] depends linearily on 50* and their spatial derivatives. Finally, 
when Sz-^, 5u'^ satisfy the linearized Hamiltonian equations of motion, including con- 
straints, we find from (15.21) and (15.41) that 

d,kl [6z,, zs] = 0, dokl [6z,, zs] - a.tgl = 0. (5.5) 

At a fixed time t = x*^, consider a closed 2 dimensional surface S, dS = 0, for instance a 
sphere with radius r and define the surface charge 1 -forms by 

^QeA^^s^Zs] = j) d^Xikl^[5Zs,Zs], (5.6) 

where cPxi = ^eijkdx^ A dx'^. The first relation of (15.51) implies that the surface charge 
1-form only depends on the homology class of the closed surface S, 

d^x^k2[6zs,Zs]= f d^x^k2[6zs,Zs]. (5.7) 

Si JS2 

Here Si — S2 = dT,, where S is a three-dimensional volume at fixed time t containing 
no sources. For instance, the surface charge 1-form does not depend on r. The second 
relation of (|5.5I) implies that it is conserved in time and so does not depend on t either, 

j^fQe^[6zs,Zs]=0. (5.8) 

The question is then whether these charge 1 -forms are integrable, see e.g. [|33l[34l[30ll for 
a discussion. 
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5.2 Linear theories 

In the case of linear theories, the latter problem does not arise and the whole analysis 
simplifies. One can replace (15.21) by 



^as^ = z^"-^ - d,kl[z], (5.9) 



Sz 



where S / 5z^ are the (spatial) Euler-Lagrange derivatives and kl [z] depends linearily both 
on the phase space variables z"^ and their spatial derivatives and on the gauge parameters. 
One then uses (15.11) directly to eliminate 9o(7q-£:°) from the time derivative of (15.91) . to get 



d,. 



dakl 



do 



+ 5eZ J^+^eU^. (5.10) 



For gauge parameters that satisfy 

5,^z^ = Q = 5,X. (5.11) 

one then arrives at 

d,K^[z\ = -7„£°, doKAA = s\^[z,u] - (5.12) 
For a solution z^, m", the surface charges 

Q,^[z,] = (l)^d''x,kl[z,l (5.13) 



are again independent of r and t. 

When this analysis is applied to the Hamiltonian formulation of Pauli-Fierz theory, 
one finds the standard expressions 

kl[z] = 2e™7r™ - - 5™9")/i^„ + h^ni^'^^d' - 5™a'")e^, (5.14) 

while the only solutions to (15.1 II) are = —^[fj.i']^'^ + o-fi, for some constants a^, = 
—uj[uij] ■ In this context of flat space, Greek indices take values from to 3 with = 
Indices fi are lowered and raised with 77^1, = diag (—1, 1, 1, 1). 



5.3 Electric and magnetic energy-momentum and angular momen- 
tum surface charges 

The previous analysis is not directly applicable to our case since we do not have Dar- 
boux coordinates and the Poisson brackets of the fundamental variables are non-local. 
In particular, the gauge transformations (13.171) do not allow for non trivial solutions to 
5e^z^ = 0. We also have to keep the sources explicitly throughout the argument, because 
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A ^ applied to localized sources will spread them out throughout space and we need to 
check that we are only dropping terms that indeed vanish outside of the sources. 

In the presence of the sources, the constraints 7;;^ = {Hi^, T^i±) determined 

= ^an^ - (167rG)Tl, = - (167rG)T°o. (5.15) 

Instead of (15.91 ). we can write 

- (167rG)(Tlr™ + T°or^), (5-16) 

where 

kl[z] = 2Cea6vr^™ - C^{5^^d' - + - 5™^™)^^. (5.17) 

Consider now gauge parameters (x) satisfying the conditions 

1 can I an cam n cam am fa± 

lis - u - Co^^ - c > (5 18) 

The general solution to conditions (|5.18l) can be written as 

= -^[H^^^ + (5.19) 
for some constants aj^, ^[^^^ = It follows in particular that the surface charges 

1 f .2. 



Qes[^s] = j^(pd'x,kl[z,], (5.20) 



s 



do not depend on the homology class of S outside of sources. 

Assuming A invertible, the equations of motion associated to Ct = j^i^^H + -^"^ 
imply in particular that 

dohl,^ = dmK + dnu^m " '2^''' ^H,^n + e^'S^^AH, 

+ (1671^)6"" (a-1 (OT,)^„ + ^A-'drnenp.d'dkT,''^ + ^A-^d^emp.dPdkT,'"^ , (5.21) 



_i(5-A-9'"9")(2n, + ^Aa). (5.22) 



By direct computation using the equations of motion, one then finds 

doK^Zs] = {lQnG){CsTn - d,}^^^[zs,us\, (5.23) 
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with 

ijk 



'uli2na + ^ACJ - ir{/\Han.k ' Omd' Hark) 



(5.24) 



where t<;^„ = tu'^'^efcmn- The surfaces charges (15.201) are thus also time-independent out- 
side of sources. 

Finally, the surface charges are gauge invariant, 

kl[S,z] = d,r^:^^l (5.25) 
= {K'd'vi + 2vidTs^ + C^d^vl - — - j)) - ^^''u^tvi- (5.26) 
Defining 



1 
2 

we get for the individual generators 



QeA^] = -u;;,Jr-a';,P/:, (5.27) 



(167rG)P,^ = - / d^x^d"'AC,= (f (fx^{dnhr-d^K), (5.28) 
J s°° Js°° 

(167rG')P; = 2/ d'x^eabAH''^"' = -2(f d^x^eabn'"'^, (5.29) 

'5° 



(167rG)jf = 2/ d^Xmeab{AH''"'^x^ - AH^"^^x^) (5.30) 

= -2 / (i'x^ eafe (7r'''"^x' - vr^'^'x'^) , (5.31) 

(167rG)J„^'= = / 6/2^;^ (ACT'"'^ - a^AC^x'^) (5.32) 
Js°° 

= I d\m [{dnhr - d'^K) x^ - /^r' + /ia5"'] • (5.33) 

The only non-vanishing surface charges of the dyon sitting at the origin are 

= Ma. (5.34) 
As expected, they measure the electric and magnetic mass of the dyon. 

For later use, we combine A;*, kg''^ into the n — 2 forms kg[z, u] through the following 
expressions in Cartesian coordinates, 

I. _ l.lfMu]jn-2 i.[Oi] _ Zi (C oc\ 

d'^'^T = - f dT'''=+^ dr"^ ('5 36^ 
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where eap-ys is completely skew- symmetric with eoi23 = 1 and the wedge product between 
the differentials is understood. Equations (15.161) and (|5.23l) can then be summarized by 

(ifc,, ^ -(167rG)T,^, T,, =T,arrf%, ciT,, = 0, (5.37) 

where closure of the n — 1-forms T^^ follows from the conservation of the sources, the 
symmetry of the energy-momentum tensor and (15.191) . 

Remark: In fact we have checked here that the standard expressions for surface 
charges in Pauli-Fierz theory, when extended in a duality invariant way, have all the ex- 
pected properties. More interesting would be to develop the theory of surface charges 
from scratch in theories where the Poisson brackets among the fundamental variables are 
not local to see if the ones we have found exhaust all possibilities. From the preceding 
discussion we see that pseudo-differential operators will play a crucial role for a dis- 
cussion of these generalized conservation laws, as they do in the discussion of ordinary 
conservation laws for evolution equations of the Korteweg-de Vries type for instance. 



5.4 Poincare transformations of surface charges 

Suppose now that z^, m° solve the equations of motions for the conserved sources T^^{x). 
Let z'^, u'" be the solution associated to new sources T'J^^{x) related to Tj^^ix) through 
a (proper) Poincare transformation, x'^ = A.^^^x" + with |A| = 1, 

Tr{x')=k\k''pTf{x). (5.38) 

For instance, starting from the conserved energy-momentum tensors (14.41) of a dyon sit- 
ting at the origin with world-line = 5qS, one can obtain in this way the conserved 
energy-momentum tensors of a dyon moving along a straight line, z'^ = u^s + a'^ with 
m'^, a'^ constant, n^w^ = — 1 and s the proper time, 

Tr{x') = M^u^ j dX6^^\x' - z'W)^ = Ma^6^'\x'^ - z"{x')). (5.39) 

when A'^o = u^- 

Assume then that the Casi^) transform like vectors 

Csix') = A'-^Csix) = -{AiUaA-'xr + (Auj^A-'b + AaJ^ (5.40) 

which implies that the T^^ are closed Poincare invariant n — 1 forms, 

T,', (x', dx') = Te^{x, dx). (5.41) 

We can then use the following variant of the tube lemma. Suppose that at fixed time a;°, 
T°'^(a;)^"^(x) has compact support and that there exists a tube, i.e., a space-time volume 
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W connecting the hypersurfaces n : K = and il' : K = x'^ = A^^x'^ + 6° with K 
a constant such that, dW = -n + T.lf nothing flows out through T, T^^ = 0, it 
follows from Stokes' theorem and (|5.41l) that 

[Tss= [ T,. = I T;,. (5.42) 

Jv. Jo.' Jo.' 

If we now compute the surface charges for a large enough sphere S at fixed x^ containing 
both r°^(x)C''(a;) and T^°^(x)^f'^(x), it finally follows from (I5l6l) or (l537l) that the 
surface charges evaluated for the new solutions z'"^ are obtained from those of the old 
solutions through 

Qe'K] = Qe[Zs]. (5.43) 
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A Poincare generators for Pauli-Fierz theory 

In this appendix, we assume that the canonical variables vanish sufficiently fast at the 
boundary so that integrations by parts can be used even if the gauge parameters do not 
vanish at the boundary. 

In the Hamiltonian formulation of general relativity f25||, the canonically conjugate 
variables are the spatial 3 metric gij and the extrinsic curvature vr*-' . The constraints are 
explicitly given by 

= (^vr™"7r^„ - ivr^^ -^R, = -2V (A.l) 

The associated generators of gauge transformation if [^] = J d^x {'H±^^ + Hi^'^) satisfy 
the so-called surface deformation algebra [l35l[36l. 

{H[^],H[r^]} = H[[^,7^]sdI (A.2) 
[^,v]sD = ediV^-V%^^, (A.3) 
[^^vfsD = 9''{^^djr^^-V^d,^^)+ed,v'-V^d,e- (A.4) 
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When the parameters /, g of gauge transformations depend on the canonical variables, 
(lAlU) is replaced by [I37l 



{i7[/],i^M} = i/[fc], 

k = [f,g]sD + Sgf - 5fg - m, 

\f\g^{x')}n^ix') + {f^,g^ix')}n,ix') 



m 



m 



d^x' 
d^x' 



{f\g^{x')}ndx') + {f\gK^'m,{x') 



where 



Dijki = '^^(dikgji + gjkgu — QijOki), 



(A.5) 
(A.6) 

(A.7) 
(A.8) 



(A.9) 
(A. 10) 



(A.ll) 



Let gij = 5ij + hij and consider the canonical change of variables from gij , tt'^' to = 
(hij, vr^'). We will expand in terms of the homogeneity in the new variables and use the 
flat metric Sij to raise and lower indices in the remainder of this appendix. Furthermore, 
Greek indices take values from to 3 with yU = (J-,«). Indices are lowered and raised 
with r]^tj = diag(— 1, 1, 1, 1) and its inverse. Let u5p 

To lowest order, i.e., when gij = Sij, the vector fields 



^p{uj,a) 



(A. 12) 



equipped with the surface deformation bracket form a representation of the Poincare al- 
gebra USD, 

[^p(cDi,ai),^p(u;2,a2)]?D = ^p([u;i, u;2], 25ia2 -'^201)- (A.13) 
For the gauge generators, we find H[^] = H^^^ [^] + if (2) ^ ^(3) _^ ^ where 

1/(1) [^] = j d^x {-2d^nijC + {d'd^hij - Ah)^^) (A. 14) 



d^x 



(A. 15) 



are the gauge generators associated to the constraints (12.31) of the Pauli-Fierz theory. Be- 
cause 
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we have to lowest non trivial order 

{H^'^[^],H('\v]}=H^H^,v]Pj,]. (A.17) 

This means that H^'^^[r]] are observables, i.e., weakly gauge invariant functional. 
One can use integrations by parts to show that H^^^ [^p] = 0. It then follows that 

{H[^p],H[r^p]} = {H^''>[UH^'^[Vp]}+Oiz'). (A.18) 

For vectors ^p{uj,a),r]p{9,b) of the form (IA.12I) . the first term on the RHS of (IA.16I) 
vanishes on account of (IA.13I) . To lowest non trivial order, (IA.2I) then implies 

{H^'^iUH^'^iVp]} = H^'\[^p,Vp]Pd]+H^'^[[^p,Vp]%]- (A.19) 

The generators H^'^^[^p] equipped with the Poisson bracket thus form a representation of 
the Poincare algebra when the constraints of the Pauli-Fierz theory are satisfied. Explic- 
itly, the term proportional to the constraints is 

H^'\[^,V]%] = -2 I d'xd^7r,,h^\^^9\ - v^u;\), (A.20) 

while 

= -2d, {iT^'^Kk) + -K'^dih.k (A.21) 

+ -^duK^d^h'^ - ^dkh^'d^hij + ^dihdjh'^ - ^dihd'h 

+ di (^hd^h - K'^d^hij - hid.h''' - h'^dih + hi^^d'Kj + hii^d'h^^^ . (A.22) 

Isolating terms proportional to the constraints, we find 

if(2)[^] = j {nmh"'% + \nhe^ + H^^\il (A.23) 
T^f ) = --K^^djhk, + duh,, - d,h,k), (A.24) 



+ di (^h'^d^hij - K'^dih + hi^^d'hij + hii^d'h^^^ , (A.25) 

with [^] = J (nPc + ^J^^^) • On account of (IXTTl) and the analog of (Ia31) 
for it follows that 



{mUH^%p]} - H^'\[^P,Vp]Pn], (A.26) 
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where ~ means an equality up to terms proportional to the constraints Tim-, Ti.± of Pauli- 
Fierz theory. Note that the functionals H^"^^ [^p] and H^'^^ [^p] generate transformations of 
the canonical variables that are equivalent because they differ at most by a gauge trans- 
formation of the Pauli-Fierz theory when restricted to the constraint surface. 

The generators for global Poincare transformations of Pauli-Fierz theory can then be 
identified as 



^fiu — ^fiu, 0'± — O-X, Cli — CLi + ^^±iX . 



Indeed, differentiating (|A.19I) with respect to b± gives 



(A.27) 



(A.28) 



When combined with (IA.19I) and (|A.27I) . this shows that, on the constraint surface, the 
generators Qg(^) o,) are conserved and satisfy the Poincare algebra. 

Finally, we can further simplify the explicit expression for H^"^^ [^p] by using linearity 
of ^p in a;* and integrations by parts to show that 



+ -dihd'h + di{hdih'^ + h^'d'Kj) 



ip- (A.29) 

The expansion of the gauge transformations (IA.9I ). (lA.l II) gives to first order: 



(0) 



Wd^h + Ah'' - dkd'y^ - dkd'h"" - 5'' {Ah - dudih^^) 



(A.30) 
(A.31) 



(A.32) 



B Riemann tensor and canonical variables 

By following [[Tl [TT]| (up to conventions), we show in this appendix that the duality ro- 
tations that we have defined coincide on-shell with the standard simultaneous duality 
rotations among the (linearized) Riemann tensor and its dual together with those for the 
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electric and magnetic conserved sources. We do this by showing how the covariant Rie- 
mann tensor is expressed in terms of the canonical variables. This gives us the appropriate 
generalization of the Gauss-Codazzi relations in the case of both electric and magnetic 
sources. 



B.l Covariant equations in the presence of magnetic sources 

Our conventions are as follows. Define eai...a„ = e'*^ to be totally skew-symmetric 

with ei...„ = 1. The Levi-Civita tensor is £ai...a„ = \/Meai...a„- Indices on this tensor are 

/No- 
raised with the metric, which implies that £"1 = i^gai -cn where a is the signature 

\/\9\ 

of the metric. Our convention for the dual is *tUai...a„_j, = ^oj^^'"^''Sb^...bj,a^...a„_p- 

In flat Minkowski spacetime, traces are taken with the flat Minkowsi metric "q^^j = 
diag(— 1, 1, 1, 1). Start with the "(linearized) Riemann tensor" Rj^^p^ = Ri_wpa, with only 
symmetry properties skew- symmetry in the first and last pairs of indices, 

RfjLvpcr Ri/fipa Rpuap- (B.l) 

Its double dual (cf. MTW (Ml) 



'^fiupcr — ~Gpvpcr ^ ^M'' ^ Raff^ pen (B.2) 

has the same symmetry properties. If ordinary duals of R^i,ptj aad-Gpypa are taken with 
respect to the last pair of indices, and we define 

1 

"pupa — J^pupa — 2" 
1 



^uupa — Rpvpcr 2^^'^ ^^O-Ppcn (B.3) 

^pvpa — ~^pi/p(7 — 2 ^aPpm 



Rpvpcr c^^pu ^aPpcTi 

'^iwnrr ~ ^/J!/p(7 ~^pv ^ Ra0p(n (B.4) 

we have 

^pupcT ~ ^ I J^jbpupa^ '^pvpa ~ ^ ~^bpupcT- 



^tiupcT ~ ^ \ ^jbpupa^ '^avpa ~ ^ ~^bpupcT- K^-^'J 



The 36 independent components of the Riemann tensor can be encoded in 
The Ricci and Einstein tensors are defined as 



pa R'"^ t^'^ _ryaa pa _ pa /tj n\ 

^pu ^ pauy ^ pu fJ'Oii' -^up D If^^ ■ \^-') 



1 

-'I 

If electric and magnetic conserved sources are T^'^ = (T^'^,Q^'^), with T^'^ = T^^ 
symmetric, dpT^^" = 0, the duality rotations are defined by 



af /\/ra r?^ rpaf /\/ra rjib 

pupa - b^pvpa^ J- pv - b-l- puJ 



M,aM\ = 6ab. (B.8) 
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For a tensor K^"', let K^"' = K''" - \r]^'''K. 

It has been shown in [[TT| that the duality invariant equations of motion are 

Gr = SttG T,^'' ^ R%^, + i?^,,, + i?;:,,, = SttG e^^s,.,. tJ^. (B.9) 

They imply in particular that, on- shell, the tensors R^^^p^ are symmetric in the exchange 
of pairs of indices and that i?^^, G^j^ are symmetric. Furthermore, the Bianchi "identities" 
read 

^ 9^/?!^^^ = SttG (s-^Tf - 9^Tf ) , (B.IO) 
while the contracted Bianchi identities are 

9,Gr = 0. (B.ll) 

Let 



Defining 



[Vt^pVua - VpcrVi^p] ■ (B.12) 



Rpupa - R^pa - T^^^^^paapK^^^LU ^ [-RfeAr] , (B.13) 



the tensor -R^yp^. is skew in the first and last pairs of indices, satisfies the cyclic identity 
because s'^"'"^ R^^^^^p^ = £^''P''\e°'^epaaf3K^'^ pJ^^^lR^xr] and, as a consequence, is also sym- 
metric in the exchange of the first and last pair of indices, -R^j^p^. = Rp^pu- The associated 
Ricci tensors R^;^ = R^^ - \t''^£uapaR'{r is then symmetric, Rl^ = R^^. It follows that 
R^^ = R?^„) and R?, = \e"-^£y„paRb°' ■ The Weyl tensors are then defined as usual in 



terms of R^^p^, 



Cpupa = Rpupa - K 1^vpAR\t]^ (B.14) 

and satisfy all standard symmetry properties: skew-symmetry in the first and last pairs 
of indices, tracelessness (because R^^ = K^'^ ^^yp^lR!^^]), the cyclic identity (because 
ei^P" K'^'^ p^^p„[R^^] = 0), which implies also symmetry in the exchange of the first and 
last pair of indices, 

^pvpa ^ upper ^puapy \IJ.1J) 

crp^ = o, e^'^'-^CpV = 0, c;^p^ = c;^p.. (b.i6) 

As usual, the 10 independent components of the Weyl tensor can be parametrized by the 
electric and magnetic components = {Emn, Bmn), symmetric and traceless tensors 
defined by 
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Putting all definitions together, the relation between the Riemann and Weyl tensors is 

R^pa = Cl^vpa + K^'' piupalRlr] + Ij^^'''' ^ puc^K^^ pu [Rbl^Xr)] (B.18) 
= pa + K^'' t^ypARlxr)] + T^^^^^paa^K^^pu ^[RbXr]- (B.19) 



In particular, it follows that the 36 independent components of the Riemman tensor Rj^^p^ 
can be parameterized by the 10 independent components of the Weyl tensor Cj^^^^, the 16 
components of the Ricci tensor R\^, and the 10 components of Rfxr)- 

If we define 

pa pa ^am ^ mjk r)a -i-ia pa i ca /tj tax 

'^mn ~ -'^0{m\0\n)y ~ ^ 0[j|0|fc]' '"^mn ~ '^(mn) ^ '~'mn 

the parameterization consisting in choosing the symmetric tensors 7^^„ (24 compo- 
nents), JF^, (6 components), and -R[^,^](= *R^^i,y\) (6 components) is more useful for our 
purpose. That all tensors can be reconstructed from these variables follows from the fact 
that 

pa r) aft -77 pa ca r)a -r>a ca /t) 

This means that the symmetric part of the Ricci tensors can be reconstructed from the 
variables. Since the antisymmetric parts belong to the variables, so can the complete 
Ricci tensors i?^^. Using now (|B.18I) and definitions (IB. 171) . (|B.20I) . (IB. 121) . we find 

E^nn = + T^ln) - ^(^" + TZ^) . (B.22) 

It follows that the Weyl tensors and then, using again (IB. 181) . the Riemann tensors can be 
reconstructed. 

In terms of the new parameterization, the equations of motion (IB. 91 ) read R[^,y] = 

and 

-2e"^.Fb„ = SttGTS^, (B.23) 
^7^" = StcGTSo, (B.24) 

T^ln - £'Ln + 5mn{£'' " ^7^'^) = SttGT^^. (B.25) 

Using these equations of motion, the Bianchi identities (IB. 101 ) are equivalent to 

1 
2 

2e"'9oJ-b™ = + e^.fc.F"'^) - d^S^ (B.27) 

(B.28) 



d^e^kmJ""^ + 7^^,) = -9,7^^ (B.26) 
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B.2 Canonical expressions 

We will now express the Riemann tensor in terms of the canonical variables in such a 
way that the covariant equations (IB.23I) - (|B.28I) coincide with the Hamiltonian equations 
deriving from (14.21) . 

From the constraints with sources, we find 

J-;^ = ^Ad-H:^^. (B.30) 

Assuming A to be invertible, which we do in the rest of this appendix, TZ"- and C"^, 
respectively JF^ and (9"iJ^^ determine each other. By taking the divergence, the Bianchi 
identity (IB.26|) implies that 

Similarily, the Bianchi identity (IB. 271) implies in particular that A£°- — d'^d'^S^^^ = 
e"-''doAd'^d'^Hbmn- When combined with (IB. 251) . the equations of motion following from 
variation with respect to C"' read 

Ia^C" + e'^'AdoiAHh - d^d^H,mn) + 2A^n'^ = AS'' - a'"a'^(7^^„ - Cn)- 
When combined with the previous relations, they imply that 



2 

= -^e^'doAiAH, - 29'"9"i^6„„) + A^n\ 

The rest of the Bianchi identities (IB. 261) . (|B.27I) are taken into account by applying a curl. 
This gives e^•^*9,9'=7^ffc = \A{Ad^El^ -d''d"'d'''El^^)w\d.e'''''d,d^E'i^ = e''''2e''^dodsJ^u- 
QTQkjra ^ ^jrar_ ygt anothcr curl gives afc(9™9"7^^„ - A9"7^^„ = ^ekird^A^d^'H^'- and 
9fc9"^9"£^„ - Ad^S^^ = 2e''''doidkd^J^bn - AJ^bk) + tkw&AJ^''^ Using the previous 
relations we then get 

= -\dk^'C'' - ^eMrd'Ad-HZ, 

d^El^ = e'^'doAi-^dkHb + d^Hbkn) + d^An'^ - ^e^rd'Ad^HZ. 

The equations of motion following from variation with respect to are then identi- 
cally satisfied. 

Defining = 7?.,5^„ — and using definition (12.121) of V^^ combined with 
(IB. 251) . the equations of motion following from variation with respect to read 



^abdo 



2 {V^'^H')^^+d^AA'^+dnAAl+^{6mnA-d^dn)C'] -eabA{dmn'^+dnnl)- 
- 2A^Hl^ + 5„„A2/7» = -e„,,a^Dl - e^,,^^!)^^. (B.31) 
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Taking into account definition (12.121) and previous relations, we can extract 



b 

n 

1 



\ I mn 2 " 



^"H^n + ^S^n^'H-- (B.32) 



In order to extract the remaining information from (|B.31I) . we first apply 5™" A — d"^d" 
to get 



and then a divergence 9"^ giving 



0, 



(B.33) 



eabdoi^Ai - enpgd^d'Hl") = e,,An^ + 2A5'=ffi - ^dnAH'^ - d^d^d'E'^i. (B.34) 



We can now inject the latter relations into (|B.31I) and use (12.141) . (|2.6I) to get 



-i6^^A + dmdn)AC\ 



(B.35) 



(B.36) 



Injecting into the second form of the last Bianchi identity (|B.28I) and using previous 
relations gives 



e^.nd'^AH)^ + e,,nd'^AH^ + ^(5,,A + did^)AC 



bn 



. (B.37) 



Identifying the terms with time derivatives gives 
1 



7^" 



Signd'AHp + e^gnd'AHr + -{S^,A + did,) AC 



al 



. (B.38) 



The terms without time derivatives in (IB.37h then cancel identically. Together with (IB. 361) 
this then finally gives 



= -eabdoin"^' - ^5'%') + d'd^n, - U^^'d^d^ d^^ip - le^^^'d^d^d^aip- (B.39) 
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B.3 Riemann tensor for linearized Taub-NUT 

Following for instance [|39l section Ai.2 and using a regularization in Fourier space, we 
find for the gravitational dyon at rest at the origin discussed in section l42l 

7^^. = GM« [^6.,6%x) + ^(5., - ^)] , (B.40) 

^- ■ = GM^ [f6.,6%x) + ^-^{6, - ^)] , (B.41) 

where 77 (r) is a regularizing function that suppresses the divergence at the origin and is 1 
away from the origin. We then find 

R^o = GM'^ATTd^ix), ^ = GM^AnS^jS^x), (B.42) 

E:^=GM^^M^6,-'-^), (B.43) 

and all other components of -R^j, vanishing. For the Riemann tensor, this implies 

^0.0. = GM^ [fs.,S^x) + ^(5,- - ^)] , (B.44) 

Rh, = -e^^.'GA^r [^6u6\x) + ^(5. - ^)] , (B.45) 

with all other components obtained through the on-shell symmetries of the Riemann ten- 
sor. This is the usual Riemann tensor for the linearized Taub-NUT solution. 
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